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The advent of large-scale computing equipment has made these concepts feasible. 

A logical extension of the technique considered here would be to put an additional loop 
around the object dynamic system. This loop would determine a model structure with 
which to estimate parameters. These parameters would permit state determination. This 
would approach the advanced objective previously mentioned. This research will be con- 
fined to those cases for which a rational model structure has been chosen. 

To perform a linear analysis, a mathematical model structure must be contrived or 
selected. Many considerations contribute to the selection of a rational model structure. 

The first step is to define a response matching criterion and then choose a model struc- 
ture capable of fulfilling the criterion. This fulfillment, generally checked by simulation, 
is accomplished by time response comparison, frequency response comparison or, more 
generally, a combination of both. The inclusion of small signal nonlinearities, such as 
stiction, windup, hysteresis, and deadband, is dictated by the dynamic effect on the 
system as ascertained by the response matching criterion. A central feature of model 
structure selection is system dimensionality truncation. Dimensionality truncations for 
model selection are of two types; those due to modeling complexity, and those to 
reduce dimensionality of an already chosen model structure. Modeling complexity is 
necessarily broken at some level since most complex system models could be made 
infinite in extent. These truncations are made on the basis of insight, feel, experience, 
and logistics of the computational equipment available. A rational' choice of model struc- 
ture may simply be the exclusion of dynamic effects in some frequency regime for which 
the control and/or the plant are nonresponsive. A more complicated’ scheme consists of 
including only coupling dynamics. The implementation of the coupling dynamics scheme 
is straightforward although sometimes computationally difficult. A subcomponent repre- 
sentation is determined first as an isolated system and then compared to the subcom- 
ponent representation in the closed-loop system. If the pole-zero representation moves 
more than some judgmental amount, that subcomponent must be included in the overall 
system dynamics. If the subcomponent’s dynamics may be discarded, the steady-state 
contribution is accounted for algebraically. A simple state variable criterion is to eliminate 
those states whose derivatives remain less than some judgmental amount. Many other 
schemes may be devised as well as combinations of these schemes. The ultimate criterion 
is the satisfactory working of the finally designed and analyzed system. 

The purpose of the research may be summarized as providing the “best” repre- 
sentation of the system linear model, or Jacobian, for a given configuration. Best repre- 
sentation means the best model attainable under a qualitative judgment involving 
accuracy, measurement inaccuracies, and system disturbance. The principal objective is 
the assessment of the validity of the mathematical model used to design a subject system. 

While there are other techniques providing the same information, the proposed 
technique reduces the system history required for solution. Most estimation techniques 
avoid the partially measured state vector because the standard approach is to adjoin the 
unmeasured states to the parameter matrix. The result is that an n-squared problem has 
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CHAPTER l 
INTRODUCTION 


Statement of the Problem and Objectives 

An advanced objective of control system theory is to build a learning device so an 
unknown system is directed to some goal by the device. First the identification of the 
system is implemented. Based on the identification, the state of the system is determined. 
An appropriate control stimulus results in a response that satisfies some rational figure of 
merit. This is presently accomplishable in only basic systems [1]: The purpose of this 
research, is to develop and demonstrate a technique, consisting of several concepts, that 
permits simultaneous calculation of the state and a reasonable facsimile of the plant. 

These concepts center upon state reconstruction and parameter estimation. 

There is a certain amount of literature using these concepts [2,3] , but only 
recently have the combined concepts of state reconstruction and parameter estimation 
been exploited [4] to provide information simultaneously of the state and the system 
representation. The literature is confined to linear autonomous systems, while this 
research will attempt to extend the developed techniques to nonautonomous and non- 
linear systems. The technique of Reference 4 requires a Liapunov function of the 
unknown system, which is possible for linear systems. The advantage is the synthesis of a 
globally convergent scheme. The disadvantages are that the Liapunov function may not 
exist for nonautonomous and nonlinear systems. The method of this research uses a steepest 
descent of gradient type method. A disadvantage of a gradient method is that initial 
estimates must be close. However, in practice, the system is reasonably well known and 
this disadvantage is not overwhelming. , 

This paper will attempt to apply the combined techniques of parameter estima- 
tion and state reconstruction to the measurements of nonlinear physical systems. In 
practice, linear systems do not exist. However, there are regions of operation on which 
any system exhibits nearly linear behavior. The limitation may exist that the region of 
linearity is too small or that expected excitations will drive the system out of its linear 
region. In any case, a comprehensive control or system analysis begins with a linearization 
of the subject system. For many applications, the analysis either forms a basis for design 
or provides a rationale for redesign or alteration of the system. A fundamental, but some- 
times unanswered, question is “how good a representation of the system is this linear 
model?” In many fields and applications, an a posteriori analysis is undertaken to assess 
the mathematical modeling accuracy. Generally, this consists of a manual Iterative assess- 
ment until some degree of accuracy is achieved. 

Parameter estimating state reconstruction can be applied, as a black box, to a 
system’s measurements to provide a real-time assessment of the current linear model of 
that system. By assessing the response from different pieces of hardware, statistics may 
be compiled as' to the spread of that system’s operation. An analysis of extreme condi- 
tions provides an assessment of the sensitivity of the system to a real environment. As a 
result, a quantitative assessment of the system design analysis is possible. 



been expanded to an n-squared plus n problem. On the other extreme, the Liapunov 
approach becomes difficult for systems of greater than single output because a Liapunov 
function must be contrived for each of the outputs which also are coupled. The proposed 
technique retains the n-squared dimensionality of totally measured systems. 

Figure 1-1 is a diagram of the system description. A state reconstructor is used to 
recover the unmeasured states. The state reconstructor is used with the reference model 
to form an error for the estimation process. Establishing validity of using the state recon- 
structor output and addressing coupled systems dynamics are the principal concerns to be 
analyzed. 


Chapter Description 

Chapter I is a delineation of the area of research with an identification of prob- 
lem areas and overall objectives. Chapter II will develop the observability theorems for 
the systems to be considered. This will be followed by developing and presenting the 
necessary observer or state reconstruction theory to support the research. The treatment 
will be subdivided to treat linear time invariant, linear time varying, and nonlinear sys- 
tems. Chapter III will address parameter estimation with state reconstruction. The prob- 
lem of model structure and model matching criteria will be analyzed. 

Chapter IV will develop a dynamic model of the CTL-V Space Shuttle POGO test 
facility to be analyzed by the technique of this paper. CTL-V is a particularly good 
example since fully half of the involved parameters are unmeasurable. A nonlinear model 
will be developed and the linearized equivalent will be analyzed dynamically at the rated 
power level operating point. Chapter V will design the necessary observers, to permit 
application of the technique, and present the results of the simulation and analysis. 
Finally the technique will be summarized considering its advantages, disadvantages, and 
unique characteristics. 
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Figure 1-1. Block diagram of the system. 
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CHAPTER It 

OBSERVABILITY AND OBSERVERS 


Paramount to any closed-loop control consideration is the measurability and/or 
observability of the object dynamic system. An unobservable system may be controllable 
only in an open-loop sense. Therefore, this research will be confined to closed-loop con- 
trol and thus observable systems. Observability will be dealt with in detail, particularly 
those aspects pertaining to linear constant coefficient, linear time varying, and nonlinear 
dynamic systems. Observability is a required condition for the state reconstruction 
process and is included for completeness. 

Many unobservable systems may be recast in a form tractable to the technniques 
of this research. The procedure consists of partitioning the system into observable and 
unobservable parts. The partitioned observable part may then be handled as an observable 
system. If the unobservable partition interacts with the observable partition, the inter- 
acting elements may be treated as disturbance inputs to the observable system. The par- 
titioning process can be accomplished by means of a transformation to controllability- 
observability canonical form [5]. 

A pedagogic examination of the relationship between controllability and observa- 
bility will aid in the development of the requirements for observability. If x is an n-space 
representation of the system, y is an m-space representation of the measurements of that 
system, and Xq is the initial state, then the controllability problem may be defined as the 

existence of a solution from xq to a desired state, Xp The observability problem is 

defined as the existence of a unique one-to-one mapping from x to y. This has been 
elucidated by Kalman [5] as the principle of duality. The principle of duality [6] 
depends on the uniqueness of the solution and the mapping. This principle applies to 
linear constant coefficient and linear time varying systems. However, for nonlinear sys- 
tems, existence is not necessarily uniqueness and the principle does not apply [7], It 
remains to develop the conditions for assuring, a unique one-to-one mapping from x to y 
for the various dynamical systems. 


Observability Theorems 
Theorem : (Linear time invariant) 

The system 


x = Ax 
y = C T x 
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where x is an n vector, y is an m vector (m < n)„ A is n x n, and C 1 is m x n, is com- 
pletely observable if and only if the composite n x mn matrix 

[C, A T C, . . . A Tl>1 C] 

is of rank n. The proof of the preceding is given in many texts [6,8] . 

Theorem: (Linear time varying) 

The system 


x = A(t)x • 
y = C T (t)x 

where the variables are as previously defined, is completely observable on the time -inter- 
val Iq < t < t j if and only if the matrix 


M(t 0 ,t j) - / ^ T (r,t 0 ) C(r) C T (r) <I>(r,t 0 ) dr 
fc 0 


is nonsingular. The matrix <h(t,tQ.) is the unique fundamental matrix satisfying 
$(t,t 0 ) = A(t) <kt,to) , = l n 


For complete observability, the above must hold for every E 0 and some finite tj > Iq. 
The proof of- thespreceding theorem- is likewise found in most modern control texts 


[ 6 , 8 ]. 


For nonlinear systems, a more precise definition of terms is required because 
existence and uniqueness are'no longer equivalent. The nonlinear system may be repre- 
sented as 
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X = f(t,x) 


f : [t 0 ,t ! ] xftCE' X E n -» E n 


( 2 - 1 ) 


with measurements 


y = h(t,x) 

h:[t Q.tjl xf2CE' x E n E m (2-2) 

The initial state x(t Q ) is in general unknown since m < n. Now assume that the £th order 

derivatives of f and h exist for every xe£2 and for every te[tQ,tj] where £m > n. Expand 
y(t) as a Taylor series 

y ( 2 ) (tO> n 

y(t) = y(t 0 ) + yCt 0 )(t - 1 0 ) + . . . + — — (t - t 0 f (2-3) 


where 


y(to> 

KtQ) 


h(x(t 0 ) 5 t 0 ) 4 h 0 (x(t 0 ), t Q ) 

3h 0 / 8h 0 \ 

(x(t 0 ),tQ) + y g" - (x(t 0 ), t 0 ) J f(x(t 0 Xt 0 ) A h,(x(t 0 ),t 0 ) 


Now define 


H(x(t 0 )) 

where 


y(t 0 ) 


h 0 ( x (t 0 ),t 0 ) 


• 

H(x(t 0 )) = 

• 

(2-4) 

_y (JM) ftp 


h£-i(x(to),to) 




The nonlinear map H(x(tg)) is called the “observability mapping” of the system. The 
system described by Eqs, (2-1) and (2-2) is said to be completely observable in £2 q on 
the time interval Etgdq] if there exists a one-to-one correspondence between the set £2q 
of initial states and the set of trajectories of the observed output y(t) for teftQ.tj]. If 
the observability map H is one-to-one £2g to H(£2q), then knowing z uniquely determines 
x(tg) so that the system is completely observable. Several publications [7,9] have investi- 
gated these conditions for global observability. 

Theorem: 

The system described by Eqs. (2-1) and (2-2) is completely observable in the set 
£2 q of initial states on the time interval [tgdq] if 

(1) 2m ” n, where n is the span of the state, m is the number of outputs, and 
£ is the 2th derivative of f and h which are assumed to exist. 


(2) The observability mapping of this system is differentiable. 

(3) There exists an e > 0 such that the absolute values of the leading principal 
minors Aj, A2, . . . , A n of the system Jacobian satisfy 


lAj t > e 


!A 2 i 

lAjl 


lA, 


1A 


> e 


n-1 


for all xeE n , then H is one-to-one from E n onto H(E n ). This result is proven in 
Reference 9. 

The development of the conditions for nonlinear observability gives visibility to a 
minor theorem that can be applied to linear systems. If the observability mapping is 
related to linear systems, it reduces to the familiar form of the condition for observabil- 
ity. Another use may be made of these results. First, consider the case where there is 
only one measurement. Now the vector z is simply 


y 

y 


" C T ’ 

c t a 

* 


« 

_y(n-l) 


_C T A n-1 _ 


(2-5) 
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which is the transpose of the observability matrix. This matrix must be of. rank n since 
the system is observable so that premultiplying by the transponse and taking the inverse, 
Eq. (2-5) may be written as 


fc T 

c t a 

T 

C T 

c t a 

-1 

C T 

c t a 


y 

y 

C T A n_1 _ 


_C^A n "*_ 


cTa 11 -^ 


y(n-l)_ 


( 2 - 6 ) 


Given the output and n-1 derivatives of the single output, the n-state vector may be 
deterministically obtained. This result is of little use because, in practice, it is difficult to 
differentiate the measurement with great insensitivity or precision. 


Observers 

Observer or state reconstructors have been repeatedly examined in the literature 
since Luenberger [10,1 1] quantified the concept. The reason for this interest was the 
advent of state variable theory [6], which organized a dynamic system in such a fashion 
that the observations of system are not necessarily the measure of the system. The sys- 
tem may be made up of n states and observed by m observations where n is not neces- 
sarily equal to m, but generally is greater than m. The observer fills a need to have n 
states from m observations.! 

The objective of the observer process is to provide reasonable approximations to 
those states that are not directly measured. Then these states are available for use in the 
implementation of a control law or strategy. Observers also find use in system estimation 
and identification. The conceptual basis for the observer lies in the process of driving an 
auxiliary dynamic system with the available outputs of the subject dynamic system. 

The state reconstructor is an auxiliary dynamic system that deterministically 
calculates the states using the difference between the real measurements and the measure- 
ments from the reconstructor. The state reconstructor is an intriguing mathematical 
phenomenon because apparently “free information” is acquired. That is, m measurements 
are sufficient to determine an n-vector state. Considerable attention has been focused on 
“reduced order observers. ” If some of the states are directly observed the system may be 
partitioned to form two related systems. These related systems consist of one measured 
and the other reconstructed. 

One of the most significant problems of observers is knowing initial conditions for 
the reconstructor system. The initial conditions of the measurements are apparent but 
these are not necessarily the states. This problem is complicated by errors in the system 
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parameters. A solution to this dilema is the “robust observer” that will be investigated 
further. These observers have the property of converging to the proper state, even though 
the differential equation coefficients are not accurately known. This is accomplished in a 
manner analogous to integrating out steady-state error. 


Observer Development 

An auxiliary dynamic system will almost always serve as an observer in that its 
state will tend to follow a linear transformation of the subject dynamic systems state. 
The design of the observer consists of incorporating that linear transformation into the 
process, thus providing an immediate and direct measure of the state. 

Let Dj be a free dynamic system describable by 


k(t) = Ax(t) 


(2-7) 


and D 2 will be the auxiliary dynamic system of the form 


z(t) = Fz(t). 


( 2 - 8 ) 


This auxiliary system will be driven by the outputs of Eq. (2-7) 


y(t) = C T x(t) 


(2-9) 


so that 


z(t) = Fz(t) + Hx(t) 


( 2 - 10 ) 


where 


H 


KC T 


( 2 - 11 ) 


in which K is a gain matrix selected to achieve some goal. Now 

z( t) - Px(t) = Fz(t) + HPx(t) - PAx(t) (2-12) 
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If 


H = PA - FP (2-13) 

then 

z(t) - Pic(t) = F(z(t) - Px(t)) (2-14) 

which has 

z(t) = Px(t) (2-15) 

as a solution, demonstrating the assertion of the preceding paragraph. Notice that Dj and 
D 2 need not have the same dimension. 

- This suggests the “identity observer” where the transformation P is the identity 
matrix. For this type of observer, Dj and Do must be the same dimension. Note that y.(t) 

becomes an estimate of x(t), F becomes A-H, and Fq. (2-10) may be rewritten as 

z(t) = x(t) = (A - KC T ) x(t) + KC T x(t) . (2-16) 

Let the error between x(t) and x(t) be defined as e(t). Now, 
e(t) = x(t) - x(t) 

= Ax(t) - Ax(t) + KC T x(t) - KC T x(t) 

= (A - KC T )(x(t) - x(t)) 

= (A - KC T ) eft) (2-17) 

which expresses the dynamics of Figure 2-1. 
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If CT and A are real matrices, then the eigenvalues of A - KC T can he made to 
correspond to the set of eigenvalues of any nxn real matrix by suitable choice of K if 

and only if (C 1 ,A) is completely observable. This has been proven in the literature several 
times, notably by Gopinath [12] . This implies that e(t) may be driven to zero arbitrarily 

fast by suitable choice of eigenvalues of the augmented system A - KC 1 . The response is 
normally dictated by a trade between accuracy and performance. If the eigenvalues are 
made extremely large negative, the system tends to act as a differentiator and is highly 
sensitive to noise and other disturbances. Some of these effects will be demonstrated with 
a subsequent example. 

The observer is easily expressed as difference equations. Equation (2-7) becomes 


x(k+l) = G(T) x(k) + H(T) U(k) (2-18) 

where G(T) = e AT , and Eq. (2-16) becomes 

x(k + 1) * G'fT) x(k) + L(T) x(k) + L'(T) U(k) 

where G'(T) = e (A-KC^)T 
T 

L(T) - j e^ A ‘ KC;T > r KC T dr 
0 


T 

/ J 

e (A-KC )r Bdr (2- 21 ) 

0 


(2-19) 


( 2 - 20 ) 


This may be represented as in Figure 2-2. If the system is not totally observable, then the 
ability to place eigenvalues is restricted. In fact, some of the errors may be unbounded. 
This does not imply a lack of system controllability but rather a lack of adequate control 
command to the state reconstructor. Only if the system is totally observable can the 
eigenvalues of the error system be arbitrarily placed. Without total observability, the 
reconstructor will be uncontrollable. 
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Let M be the nonsingular transformation that takes the system to observable 
canonical form: 


x = Mz 


( 2 - 22 ) 


Rewriting Eq. (2-16) 


Mz 

= AMz + BU + K(C T Mz-C T Mz) 

(2-23) 

Mz 

= AMz + BU 

(2-24) 

e = 

Mz - Mz 

(2-2 

c = 

Mz - Mz 

(2-26 

6 = 

(AM - KC T M) M _1 e 

(2-27) 


But 


c t m = (C r ■ O mk ) 


where C 1 is in x ink where k states are unobservable, and KC7M is n x n-k. it is 
apparent that the last k columns of AM are unaffected by choice of K. 

Nonlinear observers have been developed for several cases [13, 14,1 5J. In general, 
these observers arc highly system dependent, and are very sensitive to initial conditions 
and gains. System dependent means that the closed form observer may be developed only 
for distinct classes of nonlinear systems. Further restrictions are the conditions required 
for one-to-one mappings which assure observability. More general realizations of 
observers, characterized by Reference 13, require nonlinear gain schedules for con- 
vergence and limitations on initial conditions. 


Observers for Use in Estimation 

The principal problems in applying observers to estimation are isolating the 
dynamics of the system from the dynamics of the observer and knowing the matrix A. 
The observer existence is based on some knowledge of the matrix A. While in many cases 
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the accuracy of the reconstruction process is independent for small errors of the precision 
with which A is known, there are cases where the observer will diverge [16], Due to 
manufacturing tolerances, material acquisition, sensor accuracy, modeling truncation, and 
a myriad of other reasons, the A matrix will never be precisely known. The concept of 
ah observer is analogous' to a pole-zero canceling compensation so that sensitivity is an 
inherent design problem. 

Observers of rank less than n are known as reduced order observers. The unrecon- 
structed states are obtained directly from measurements and the system is partitioned to 
separate the directly measured states from the rest. The remaining states are recovered by 
an observer of order n minus the number of measured states. The minimal observer 
results when all the measurements are used to identify specific states. However, the 
minimal observer results in a totally open-loop observer for the unmeasured states. Philo- 
sophically, reduced order observers are attractive due to the reduced dimension. In prac- 
tice the reduced order observer only simplifies the observer design. This simplification is 
easily outweighed by certain advantages of the identity observer. A minor theorem [17] 
shows that any identity observer is robust. Battacharyya [17] defines a robust observer 
as a closed loop system, closed on the error between plant and estimate, and one that 
possesses redundancy. Measurement redundancy means that, implicitly or explicitly, at 
least one linear combination of states is contained in the measurements. 

These recent studies [16,17] have been directed to the sensitivity of observers. 
While the emphasis of these studies has been on reduced order observers, the sensitivity 
results generally apply and will be used as justification for assertions and assumptions of 
this research. The primary assertion is, if a robust observer is designed, the estimated 
states will converge asymtotically to the states even if errors exist in the estimate of the 
plant parameters. The use of identity observers removes the concern that the observer is 
not robust. Further the observers will be designed so the augmented system eigenvalues 
are critically damped. This stipulation is of little value for linear systems not disturbed by 
random inputs. However, for weakly nonlinear and perturbed systems, intuition and 
experience indicate the critically to highly overdamped roots will behave in a superior 
manner. 
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Figure 2-2. Difference equation form of the identity observer. 













CHAPTER 111 

PARAMETER ESTIMATING STATE RECONSTRUCTION 


Complex multiloop system analysis and control design are most generally pre- 
dicted on a linear analysis. Sufficient mathematical modeling is developed to assure a 
reasonable facsimile of the physical process, at least in the region of interest. Therefore, 
it is highly desirable to quantify differences between behavior, in the small, of the 
process and its analysis linear model. This quantification provides a final step in the 
design cycle, and determines expected variations from the process description used in the 
system design. These variations may he sufficiently large to dictate redesign of the con- 
trols or of the process itself. Unfortunately, all states are generally not available through 
measurement, so that the process of determining the best linear model is coupled with 
determining the states on a time history basis. The basic assumptions of this chapter are: 
that a model structure has been determined, that the system inputs are known, that there 
is no input disturbance or measurement noise, and that measurement of the state is 
incomplete. 


Estimation of Partially Measured Systems 

There are two alternatives available to solve the problem of estimating the param- 
eters of a system whose states are not all available through measurement. The first is to 
augment the parameter estimation problem with the unmeasured states. The second is to 
develop in-line estimation algorithms which permit the simultaneous calculation of the 
desired parameters and states. The former results in an increase in dimensionality of the 
problem. In general, the parameter estimation problem is an n-squared problem. Augmen- 
tation can raise this to as high as n-squared plus n. This causes a long data stream to be 
required in the claculations and may result in convergence problems due to the age of 
oldest data. The latter approach is likewise potentially an n-squared plus n problem, but , 
the data stream is the same as for a fully measured estimation problem. There are several 
methods available for the calculation of the states. The most straightforward approach is 
the state reconstructor or observer. 

Gates [13] points out that the use of standard linear observers causes problems 
due to the coupling dynamics between the state reconstruction process and the plant 
estimation process. These problems are of particular concern if the dynamics of the plant 
are fast. This robs the designer of flexibility in locating system poles and zeroes to better 
solve the problem. Some of the difficulties of observers under variations have been 
addressed by Battacharyya [16]. However, if the system is stable the robust observer 
may be designed. Robust observers are stable and converge to the actual states even if the 
system differential equations are not accurately known. The reason for this convergence 
to the actual state is quite analogous to an integrator in the feedback of a control sys- 
tem. The steady-state error is driven to zero. 
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If the system dynamics are reasonably fast and if the plant is stable, the design 
approach is straightforward. The approach may be described as a periodic sequential 
estimation process. This is implemented as an estimation update periodically occurring 
after the settling time of the state reconstructor has expired. This process may be viewed 
as a sequence of state reconstruction problems with the error in parameter estimate form- 
ing the initial condition guess for the subsequent state reconstruction. The sequence may 
be designed to provide a set of essentially uncoupled state reconstruction problems, each 
with a successively better estimate of the differential equation coefficients. .For example, 
if the system has a 0.1 sec settling time, an observer with a settling time of 1 sec is 
uncoupled from the system and generally four or five iterations provide reasonable con- 
vergence. Therefore in this example, the estimation process may be completed in 4 to 5 
sec. There is inherent design flexibility to adjust these settling times (of observer and 
estimator) for a desirable response. 

This approach appears moic desirable than the nonlinear observer-like system 
proposed by Gates. The design of Gates’ observer is much more system dependent than a 
system tuning approach. This is undesirable, because sufficient information concerning the 
system may not be available. A further, and more serious, disadvantage of Gatos’ observer 
is that a stable configuration may not in fact exist. 

The problem is then to keep the n-squared estimation problem while developing 
simultaneously the system states. Consider a robust observer compared with a best system 
estimate or system model. If the systems aie adjusted so that the overall computation 
cycle is longer than the augmented system settling time, then the stales Ijoin the first 
may be used to reestimate the dynamics of the second. The estimation problem is now 
identical to that of the literature 1 8,13,18) since the accurate states are now available for 
the estimation process 

The observer system is released with its initial conditions set to the best guess of 
the states. After one settling time, accurate states will be available for estimation of the 
dynamics of the system estimate using the method of Gates (.131 . The system estimate is 
needed to use as a trial horse against the observed system behavior. These systems may 
be successively stepped to the desired convergence or to follow a slowly time-varying 
system or a weakly nonlinear system. 

Time-varying, linear systems become more interesting and more complex. Now the 
dynamics of the estimator must be rapid enough to track the system under expected 
variations to fulfill the response matching criterion. The state reconstruction dynamics are 
chosen to be faster than the estimator within considerations of physical restrictions and 
overall stability. For nonlinear systems, one may consider first those that possess only 
small signal linearities or those that are weakly nonlinear. A well-known artifice of con- 
trol analysis is to treat weakly nonlinear systems as time-varying linear systems. This 
approximation is not without risk, since an unstable system may appear to be stable in 
this type of analysis and vice versa. If the time-varying approximation is permitted, the 
technique may be applied as with the time-varying case with the additional restriction of 
defining a region for which the response remains sufficiently weakly nonlinear. The 
example of Chapters IV and V is of this typo. 
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Stability considerations for these processes are concerned with three distinct ele- 
ments. The stability of the plant estimate is of primary importance. Parameter variations 
in the estimation process may move marginally stable eigenvalues into the right half 
plane. If this occurs special tuning techniques are required to achieve the desired accuracy 
due to the observer [16] behavior. The second element is the augmented system eigen- 
values of A + KC T . In practice this is seldom a difficult task because the elements of K 
are quite large in the interest of rapid system response. The variations in parameter esti- 
mates are normally small, in relation to the elements of K; thus, augmented system 
eigenvalues tend to be insensitive. The last consideration is the stability of the estimation 
algorithm. Many techniques exist [18] to assure the stability of this process. 


System Diagrams and Equation Development 

The equations and block diagrams will now be developed to give form to the 
method. Since most implementation schemes are digital in nature, the equations and 
algorithm will be developed in discrete form. 

The implementation has, as its objectives, the recovery of the unmeasured states 
and the estimation of the system parameters. The unmeasured and measured states are 
direct outputs of the observer previously described. These states are used directly in the 
estimation process developed by Gates [13]. Figure 3-1 is a block diagram of the system 
estimation process. The plant dynamics are differenced with the dynamics of the refer- 
ence model to form an error which is* used to determine the difference between the refer- 
ence model and the plant. The reference model is updated with the calculated difference, 
scaled by an appropriate gain. This gain is chosen for the stability of the estimation 
process. The process is repeated until the desired convergence is achieved. 

Figure 3-2 portrays the combined state reconstruction and parameter estimation 
process. The plant, the observer, and the reference model are driven by the input U. The 
observer is also driven by the measured output of the plant. An error is formed from the 
difference of the observer and the reference model. Due to the behavior of the robust 
observer, the plant estimate x will follow the plant x. The observer and the reference 
model use the same estimate of the plant. The previously determined error will be used to 
calculate the difference between the plant and the reference model. The algorithm per- 
forms this function. The reference model and the observer are updated by use of the 
calculated difference and again scaled by an appropriate gain chosen for estimation stabil- 
ity. The process is repeated until a desired degree of convergence is achieved. 

The equations will now be developed considering known external inputs. The 
plant is assumed governed by 

xCk + 1) = G(T) x(k) + H(T) U(k) (3-1) 
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where 


G(T) = e AT 


(3-2) 


and 


(k+l) T 

H(T) = f e Ar B dr (3-3) 

kT 


with measurements 

y(k) = C T x(k) . (3-4) 

A is the system Jacobian, U is the known input and B is the control distribution matrix. 
The observer is described by 

£(k+l) = G'j(T) x(k) + HjCDlKk) + L'j(T) y(k) (3-5) 

where 

G j (T) = e ( D+KcT ) T ' (3-6) 

(k+l) T 

L(T) = f e (D+KC T )r g dr (3-7) 

kT 


and 


(k+l) T 

L\ = / e( D+KcT > r Kdr (3-8) 

' kT 

D. is the estimate of the system Jacobian, and K is the observer gain matrix. 
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The reference model is expressed as 


xjdc+l) = Gj(T) xj(k) + Hj(T) U(k) 


(3-9) 


where 


G|(T) = e DT 
and 

(k+l ) T 

Hj(T) = J e Dr B dr 
kT 


(3-10) 


(3-11) 


Now, writing state in terms of the system estimate 


x(k + l) = G(T)x(k) + SG(T)x(k) + H-,(T)U(k) + SH(T) U(k) 


.(3-12) 


The error between the estimate and the state may be expressed as 


e(k + 1) = Gj(T)e(k) + «G(T) x(k) + 


SH(T) U(k) 


(3-13) 


The following matrices are formed 


xO = (x(0), . . 

. , x(-n - p + 1)) 

U0 = (U(0), . 

. . , U(-n - p + 1)) 

eO = (e(0), . . 

. , e(-n - p + 1)) 

cl = (e(l), .. 

. , c(-n - p + 2)) 


(3-14) 


where each column is the vector associated with the enumerated time point. Eq. (3-13) 
can be expressed as a matrix equation that has 
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(SG(T):6H(k)) = (el - G^T) eO) 


-1 


(3-15) 


*0 

• • 

u 0 


as a solution if the inverse of the augmented state and control matrix exists. 
The system description is now updated as 


GjCT) = Gj(T) + 6G(T) 
Hj(D = HjCT) + «H(T) 

G | (T) = G'j(T) + 5G(T) 


(3-16) 

(3-17) 

(3-18) 


Eq. (3-18) is not exactly correct but, due to the magnitude of elements of K, yields 
satisfactory results. Several alternatives exist to precisely calculate Gj(T). One method is 

to recover 6D from 5G(T) and recalculate Gj(T) from Eq. (3-6). Another method is to 

calculate as a continuous system using the numerical integration to discretize the system. 
The calculation, in this case, yields 6D. A variation of the last method would be to use 
the 6D to calculate the discrete parameters of Eqs. (3-6), (3-7), (3-8), (3-10), and 
(3-11). Precise knowledge of G'^(T) is not required due to the nature of the robust 

observer. A topic for future study is the development of a recursive solution to the pre- 
viously described estimation. 


A Second Order Example 

To demonstrate the procedure the system 




(3-19) 
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with measurement 


y = (0 



Discretizing this system results in difference equations of 


xj(k+ ly 

x 2 (k + 1), 


,-T _ a -2T „-T . 0 -2T 


2e“ 1 - e’ 


e * - e 


-2e‘ T - + 2e"^ -e"*' + 2e -22 " 


the form 



(3-20) 


(3-21) 


where T is the sampling period. 

First the observer dynamics will be developed. The error between the state and 

T 

the state estimate is dynamically determined by the eigenvalues of A + KC , or 


r 

- 


“ 



“ 

" 

0 

1 

+ 

0 

Kl 

— 

0 

Kj + 1 

-2 

-3 


0 

K 2 


-2 

K 2 -3_ 


which has a characteristic equation of 

A 2 +. (3 - K 2 )A + (2K } + 2) = 0 


(3-22) 


(3-23) 


Hence the values for 


A = 


K 2 - 3 ± y/l + K 2 2 - 6K 2 - 8Kj‘ 


(3-24) 
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are the robust observer eigenvalues and may be placed almost arbitrarily. Roots may both 
be real as K 2 = -2, Kq = 2 which has roots at -2 and -3. Figure 3-3 shows the response of 

the unmeasured state and state estimate with a perfect model. The settling time is almost 
3 sec. Figure 3-4 considers the response for a reasonable estimate of the state. Now the 
settling time is almost 4 sec. These dynamics arc too slow for application of the estima- 
tion technique. If Kj = 24, K 2 - -12 which has roots at -5 and -10 then the dynamics 

are much more suitable to the application of the estimation process. Figure 3-5 presents 
the response for these observer dynamics and a perfect model. Notice the settling time is 
of the order of 0.75 sec. Figure 3-6 shows the response for a reasonable estimate of the 
system dynamics. Again the settling time is less than 1 sec. 

Figure 3-7 shows the remarkable property of robust observers to converge to the 
state even though there are errors in the plant estimate. The settling time is on the order 
of 3 sec, and this can be improved by adjusting the observer. The observers are all 
processing data from the plant 


x = 




(3-25) 


while the reasonable estimate was 


A 

X 


/-0.2 

0.9 

\-2.2 

-3.1 


x - KC T (x - x) 


(3-26) 


The gross model error is characterized as 


A 

x 



A 

X 


- 6 . 


KC T (x - x) 


(3-27) 


which is considerably at variance with Eq. (3-25). 

The parameter estimation is delayed until reasonable convergence of the observer 
is achieved. The estimation is then initiated with the error between the trial system and 
the observer defined to be zero. Figure 3-8 presents a plot of time versus parameter esti- 
mate showing the estimation convergence. 
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Figure 3-1. Block diagram of the estimation process. 



Figure 3-2. Block diagram of the combined reconstruction 
and estimation. 
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Figure 3-7. Observer with gross model error. 
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Figure 3-8. Observer (24, -12) convergence of D to the A matrix. 
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CHAPTER IV 

CTL-V TESTING ANALYSIS 


All large space vehicles possess a longitudinal dynamic coupling of structure and 
propulsion predictably called “POGO.” The most benign stage to date is the Saturn IB 
stage which has nine tanks in a bundle and eight engines providing a maximum of statis- 
tical interaction, which results in overall system damping. The Space Shuttle, due to its 
large single LGX feedline, is expected to be susceptible to POGO. Due to this concern, an 
accumulator Is being designed for the Space Shuttle main engine as a decoupling and 
suppressive device. A primary qualification test for this device is the CTL-V Test series at 
the Rockwell International Rocketdyne Division, Santa Suzanna Test Facility. These tests 
will provide assurance as to the dynamic representation of the low pressure oxidizer 
pump and the effect and effectivity of the accumulator. 

The dynamic head rise characteristics of the low pressure oxidizer pump are non- 
linear and not precisely known. They are modeled as nonlinear differential equations 
whose coefficients are empirically determined parameters. The accumulator characteristics 
are also ill defined because of the difficulty of obtaining good test results and isolation of 
the higher frequency effects of the accumulator. The advantages of applying the pre- 
viously developed techniques are that by matching the time response the frequency 
response is likewise adjusted. That is, the linearization of the appropriate time response 
provides a frequency domain representation of the dynamic phenomena! This should be a 
“best” linear representation at that condition because the estimate is being forced to 
behave in a fashion similar to the actual system. The method is developed and modified 
to provide a neighborhood of operation of the low pressure oxidizer pump and of the 
accumulator. 

The pump modeling is defined by the Rocketdyne publication RL00001 [19} 
defining the Space Shuttle main engine-engine balance and dynamic model. The facility is 
modeled in a similar fashion. Since the lines and pumps will be essentially chilled to a 
steady state during a given test, the assumption of incompressibility and thermal steady 
state is valid. Tests have shown that there is energy trade between temperature and pres- 
sure but that these are small effects. The specific objectives of the analysis will be to 
better define the head rise dynamics of the pump and the parameters of the accumulator. 
The- accumulator parameters are characterized by electromechanical analogy. These 
parameters consist of a compliance, an incrtancc, and an equivalent resistance. The reason 
for an equivalent resistance will be apparent in the equation development. The results of 
these tests and analyses will be used in the overall Space Shuttle POGO stability analysis 
to better define system stability before first flight. 


CTL-V Equation Development 

The system to be tested in CTL-V is that of Figure 4-1. The pump speed Sqj 

may be assumed constant since the pump is being driven by an extremely high inertia 
electric motor. The constant pump speed allows analysis of the basic head rise charac- 
teristics of the pump uncoupled from available drive torque and torque required which 
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couple back into flow and pump speed. Pj is the pressure at the feedline inlet. The flow 
in the feedline DWp F may be represented as 


t 

DW FL = l” / [ (P T' P OS) ' R L DW Fl] dr 
0 


(4-1) 


The “bubble” on the pump has pressure Pqs- 


t 

P OS = r - / (DW FL - DW 0S ) dr 

l B 

0 


(4-2) 


These two elements combine to simulate the 2.5 Hz first resonance of the ozidizer feed- 
line. All damping arises from the resistive term of the feedline flow Eq. (4-2). 

Next is the low pressure oxidizer pump. The pump is assumed to have a dynamic 
gain of one. Mass continuity dictates that the pump be gain one at zero frequency. The 
pump is characterized simply as a nonlinear head rise device. The pump disciiargc pres- 
sure is 


P ODl “ P OS + H ' 


(4-3) 


The head rise H\vill be defined by use of a dimensionless parameter 


A 1 

^OPl - O— DW OS 
*0 1 


(4-3 a) 


The head rise itself is given by 

h' = BjSqj rp 0P1 (<i> 0P1 ) 


(4-3 b) 
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where PpoPl^OPl) determined from the empirical curve of Figure 4-2. The flow 

below the pump is the same as the flow existing in the bubble and entering the pump 
and is defined as 


t 

DW OS = % / [ (P ODl ' P 012> ' R D< DW OS> 2 ] dr * < 4 ' 4 ) 

1 0 

Ljq is the inertance of the fluid in the duct and the pump. R D is a lumped resistance 

coefficient combining effects of duct and pump. The pressure upstream of the accumula- 
tor is dependent on the compliance, Cq, of the duct itself, and has the form 

t 

p OI2 - / (DW 0S - DW 01 ,j - DW a ) dr . (4-5) 

0 


The accumulator is modeled analogously with a pressure change through a com- 
pliance and a flow change due to a resistance and delta pressure. The inertance in this 
case is the mass of that fluid trapped in the standpipe leading to the accumulator. The 
compliance is a lumped compliance consisting of flexure of the housing and the com- 
pressibility of the gas in the accumulator. The gas to be used on the Space Shuttle is 
GOX in contrast to helium principally used in the past. The GOX is supplied from the 
tank pressurization heat exchanger and maintains a constant level in the accumulator by 
use of an overflow port from which GOX is vented back into the feedlinc above the low 
pressure pump. That flow is not considered in this analysis. The gas-liquid interface is 
maintained by four layers of 3/8 in. Teflon balls which provide internal slosh suppression 
and help prevent the gas bubble from collapsing in the liquid. The pressure in the 
accumulator is modeled as 

t 

P A = F- / DW A dT ■ (4-6) 

A 0 


while the accumulator flow can be represented as 


30 



t 

dw A = j" / [ <P 0I2 - P A> - R A( DW A> 2 ] ^ 0-7) 

0 


where C^, L^, and R^ are as previously defined. 

PqP 2 is, in effect, an output of the system and, because of the orifice, will 

remain constant once the test conditions arc established. With one exception, the equa- 
tion definition is now complete. The exception is that to represent Pqi 2> the flow down- 
stream of the accumulator must be modeled. This is a small piece of fluid and results in a 
high frequency root. Ljq represents the inertance of that element and Rjyj is its resistance. 

The flow in this section can then be expressed as 

t 

DW OP2 = l n I ' P OP2> ' R N (DW OP2) 2 J dr • 

0 


This completes the equation development required to analyze the test. 

These equations may be rewritten as a set of nonlinear differential equations: 


R. 


nw FL = . ■ p t 


os 


DWp 


FL 


(4-9) 


OS1 


U ■ C 


DW 


B 


OS 


(4-10) 
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9r POPl 1 , „ 
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(4-11) 
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(4-13) 


Fa = 57 dw a 

'-A 


(4-14) 
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012 
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(4-15) 


R 


DW, 
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— DW, 

l n 


0P2 


(4-16) 


Linearized Analysis of CTL-V 

These equations are in turn linearized to obtain the following set of linear differential 
equations: 


1 1 R L 

AF FL = r~ AP T ' r AP OS ' 2 7 DW FL af fl 

l l 


(4-17) 


The term 2RpDWp L is normally thought of as being an equivalent damping resistance for 
an element. 


AP 0S = FT- af FL ‘ cl AF 0S (4 ' 18) 

AP 0D1 = pr af FL ■ D 1 AF 0S + D 2 AP ODl ' D 2 AP OI2 
C B 
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where 


1 


°1 = Gg ' 2A 1 B 1 S 01 ^ 


ar PO pi Rd 


opi l d 


dw, 


os 


ar POPl 1 

°2 “ A 1 B 1 S 01 lT 


OPI 


Completing the equations, 


1 


1 


AF os “ r~ ap odi - r~ AP oi2 
l d 


2R d dw os 


AF, 


-D 


OS 


AP. 


012 


AF OS ‘ AF 


'D 


-D 


OP2 


AF. 


"D 


ap a = c 1 - af a 

A 


2R. 


AF a = — AP 


012 


AP. 


dw a af a 


af op 2 " rrv AP oi2 


2R 


N 


-N 


“ DW OP2 AF OP 2 


(4-20) 

(4-21) 

(4-22) 

(4-23) 

(4-24) 


Notice that PoP2 ^ as dropped out of the linear representation because it is assumed 

constant due to the orifice. At the beginning of the chapter, the equivalent resistance of 
the accumulator was discussed. The resistance is equivalent because, to a First-order 
approximation, it is zero. In Eq. (4-23), the linear term of resistance is 2R A /L A DW A . In 

the steady state DW A is zero. Therefore, to first order the resistance effects of the 

accumulator are indeed zero. However, the resistive effects are important to the analysis 
so a small nonzero term will be forced in the analysis to assess its effect. 


33 



These equations can now be expressed in state variable form as 
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In practice, flow is an unmeasurable quantity due to instrumentation difficulties. For this 
reason the pressures are all that will be. measured. The advantages, of the scheme pre- . 
sented over most identification schemes are now obvious. Fully half of the state vector is 
not available for measurement and will be recovered with an observer designed in a 
fashion described previously. 

The measurement vector now becomes 



(4-26) 


These equations will be discretized using Fourth Order Runge-Kutta Integration, with 
iteration time sufficiently fast to assure reasonable accuracy. This provides satisfactory 
precision without cumbersome implementation. 

The test series will be operated in two ways, with and without the accumulator. 
The equations have been arranged to allow partitioning in this manner. Using data that 
reflect the rated power level test, the eigenvalues of the system without accumulator are 
as follows: 


Real 

Imaginary 

-32.45 

-224.4 

-32.45 

224.4 

-59.1 

71.7 

-59.1 

-71.7 

-5.3 

0 

-0.01 

0 
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Notice that the line frequency has dropped from 2.5 to 1.2 Hz. The pump and 
duct combine with critically damped roots, and the orifice segment has a 36 Hz 
resonance. 

Adding the accumulator, the eigenvalues become: 


Real 

Imaginary 

-81.8 

-423.1 

-81.8 

423.1 

-72.7 

77.9 

-72.7 

-77.9 

-37.5 

48.7 

-37.5 

-48.7 

-5.35 

0 

-0.01 

0 


The line frequency is essentially unchanged as is the pump and upper duct. The major 
change is in the lower duct where the duct couples with the accumulator giving the 15 
Hz resonance and the accumulator couples with the flow to the orifice which has a 69 
Hz resonance. 

It remains to locate the eigenvalues of the observer and implement the identifica- 
tion technique. 
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CHAPTER V 
CTL-V RESULTS 


The problem of higher-order observer design, analysis, and performance will be 
discussed first, then the actual design will be developed. Results will include those of the 
CTL-V facility with and without the POGO suppressive accumulator. The principal fea- 
ture is the demonstration that linear observers can be used in a state-parameter estimation 
process if sufficient care is taken in the design to insure plant-observer dynamic decou- 
pling. Techniques for desensitizing observer design and the application of parameter esti- 
mating state reconstruction will be examined. Sample rate selection and numerical 
difficulties will be addressed. Finally the curious phenomena of multiple equilibriae for 
fluidic systems of the CTL-V type will be analyzed. 


CTL-V Observer Design 

The design of observers for higher-order systems is a topic in itself. The simplest 
problem is the single input system. If the subject system is observable, there are available 
n times m parameters, of the gain matrix K, to place the n eigenvalues of the augmented 
observer system. If there are multiple inputs, the system may be recast as a set of single 
input systems and treated individually as single input systems. However, while systems 
designed in this fashion have the desired eigenvalues, the dynamics of the augmented 
system can be most undesirable because of the location of the system zeros. Undesirable 
energy trade takes place between the various component single-input systems. The addi- 
tional degrees of freedom, in the matrix K, may be used to achieve a more desirable 
overall dynamic response. The term “better dynamic response” must now be quantified. 
For the purpose of this research, better dynamic response means critically damped with a 
reasonably fast settling time, while in general the term is dependent on the application 
and the desires of the designer. 

A critically damped response is desired to eliminate or reduce coupling between 
the plant and the observer. A further precaution is to design the augmented system eigen- 
values sufficiently larger than those of the plant. This permits rapid reconstruction of the 
states and reduces the propensity of the observer dynamics to couple with the plant 
dynamics. The requirement to critically damp the observer eigenvalues means that some 
of the desired analysis flexibility has been lost and that the technique is becoming more 
system dependent. The critical damping also affects the settling time which is another 
design parameter. The settling time determines how often new estimates of the plant 
parameters can be determined. Due to the nature of observers, only the directly measured 
states are known before the reconstruction process. The unmeasured states become avail- 
able only after the system settling time has passed. This time can also be affected by the 
size of the error in the initial estimates of the unmeasured states. A settling time is 
required after each recalculation of the system parameters because the new estimate 
represents a system discontinuity when it is used in the reconstruction process. 
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A consideration that would be of low interest to most applications, but is of 
secondary interest in this application, is the augmented system eigenvalues sensitivity to 
parameter variations. This interest is of two parts: first the sensitivity of the observer 
dynamically to deviations in the plant estimate and, second, the stability of the estimated 
systems eigenvalues- due to estimation errors. Due to the nature of the robust observer, 
the dynamic behavior may be degraded as variations become large. But the system will 
converge for very large variations in the parameter estimates, as was demonstrated in 
Chapter III. In fact, the whole concept of parameter estimating state reconstruction is 
based on that property. But gross excursions can cause divergence of the reconstructor 
from the plant, and some designs are more or less sensitive to parameter variations. A 
design procedure then is to verify a low sensitivity to parameter variations. 

The second part of the problem is the closed-loop stability of the estimation 
system. The augmented system may be stable but may perform inadequately for the 
purposes of this research. If the parameter recalculation causes the observer plant to have 
unstable eigenvalues, then quite obviously the system will have inappropriate dynamics. 
Therefore, another design procedure is to analyze the sensitivity of the plant eigenvalues 
to parameter variations and the overall degree of stability of the object plant. A mar- 
ginally stable or unstable system is undesirable for analysis by this technique because of 
observer problems addressed in Chapter. II, and the overall system sensitivity. The 
observer problem is that if the plant is unstable, then the plant must be precisely known 
and represented in the observer to achieve observer convergence. However, an artifice is 
available to handle these kinds of difficulties, that is, to synthesize a feedback control 
that stabilizes or desensitizes the objectionable eigenvalues. This is a straightforward 
classical technique that provides a system with characteristics that permit analysis by the 
technique of this research, the control being included in the model structure of the new 
observer of the altered system. 

The eigenvalue placement for observers, as has been noted, is generally overdeter- 
mined. There are a variety of ways to choose the elements of the gain matrix K. The 
particular approach for this application was selected because of the manner in which the 

elements of KC^ enter the augmented system matrix A + KC^. By choosing K in the 
form 
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for the sixth order case, and 


KC 


T ^ 


K. 


21 


0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

*32 

k 42 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

K 53 


K 

0 

0 


63 


0 

0 

0 

0 

0 

0 

0 

0 



(5-2) 


for the eight order case, the augmented system eigenvalues can be- less interactively 
chosen. The system is naturally partitioned to encourage this type of gain selection. 

There is some interaction, but eigenvalue selection is more independent than if some 
more-coupled scheme were used. Values for K and the associated augmented system 
eigenvalues are shown in Table 5-1. Sensitivity results for sixth order observers are shown 
in Table 5-2 and for eighth order observers in Table 5-3. The numbers of the K matrix 

are in units commensurate with the elements of the D or A matrix so that A + KC 1 has 
meaning. If the selected elements of K are ex tremcly large, errors in the parameter esti- 
mates have little effect on the observer dynamics; however, due to the high gain, the 
observer system becomes very sensitive to noise. If the elements of K are small the 
observer becomes more sensitive to parameter estimate errors and the observer response 
becomes sluggish. These considerations enter the observer system design process. 


Physical Interpretation of the Model 

The desired output is not simply the linear model, or the matrix D in the calcula- 
tions. D must be interpreted to deduce the parameters of interest, namely the com- 
pliances, incrtances, and resistances of the CTL-V facility but, most importantly, the 
slope of the pump curve. The compliance of the bubble on the pump may be determined 
as 


C 


B 



(5-3) 
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while the feedline inertance is 



TABLE 5-1. SIXTH AND EIGHTH ORDER OBSERVER 



Sixth Order 

Eighth Order 


Value 

Eigenvalues 3 
(All Real) 

Value 

Eigenvalues 3 
(All Real)- 

*11 

1200.0 

-100.0 

1400.0 

-96.3 

K 21 

1998.0 

-219.3 

1998.0 

-138.3 

K 32 

1477.0 

-354.3 

1277.0 

-272.3 

k 42 

0.0 

-858.3 

0.0 

-544.3 

K 53 

1500.0 

-1328.1 

1500.0 

-859.7 

CO 

VO 

i4 

-1366.0 

-1496.4 

-1366.0 

-1127.6 

K 74 

— 

— 

1000.0 

-1212.2 

K 84 

— 

— 

1500.0 

-1305.8 


a. Eigenvalues have no order relationship to values. 


The line resistance is a little more difficult to recover but may be determined as 



- d 22 L L 

2W l 


(5-5) 


Next the inertance of the duct may be calculated as 



1 

D 43 


(5-6) 
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TABLE 5-2. SIXTH ORDER SENSITIVITY 



Value 

Eigenvalue 3 

Value 

Eigenvalue 3 

Value 

Eigenvalue 3 

Real 

Imaginary 

Real 

Imaginary 

Real 

Imaginary 

*11 

1600.0 

-97.2 

0.0 

1600.0 

-168.4 

87.3 

1600.0 

-184.9 

0.0 

*21 

1998.0 

-207.6 

0.0 

2498.0 

-168.4 

-87.3 

2498.0 

-185.1 

59.1 

k 32 

1477.0 

-238.8 

0.0 

1777.0 

-198.0 

0.0 

2477.0 

-185.1 

-59.1 

k 42 

0.0 

-1325.9 

0.0 

0.0 

-1337.1 

0.0 

623.0 

-1339.1 

5.6 

k 53 

1500.0 

-1378.9 

0.0 

1500.0 

-1370.0 

0.0 

1500.0 

-1339.1 

-5.6 

*63 

-1366.0 

-1508.2 

0.0 

-1366.0 

-1814.7 

0.0 

-1366.0 

-2523.1 

0.0 


a. Eigenvalues have no order relationship to values. 
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TABLE 5-3. EIGHTH ORDER SENSITIVITY 



Eigenvalues 3 

Eigenvalues 3 

Eigenvalues 3 


Value 

Real 

Imaginary 

Value 

Real 

Imaginary 

Value 

Real 

Imaginary 

K 11 

1000.0 

-98.6 

0.0 

1000.0 

-94.6 

0.0 

1000.0 

-95.6 

0 ; 0 

k 21 

1998.0 

-139.5 

0.0 

1998.0 

-146.0 

0.0 

1998.0 

-159.3 

0:0 

K 32 

1277.0 

-540.6 

0.0 

727.0 

-493.3 

-236.4 

727.0 

-492.8 

-236.9 

K 42 

0.0 

-506.6 

-212.4 

0.0 

-493.3 

236.4 

0.0 

-492.8 

236.9 

k 53 

1500.0 

-506.6 

212.4 

1000.0 

-559.8 

-393.0 

1000.0 

-555.5 

-443.8 

K 63 

-1366.0 

-860.4 

0.0 

-1266.0 

-559.8 

393.0 

-1666.0 

-555.5 

443.8 

k 74 

1000.0 

-1209.9 

0.0 

875.0 

-778.8 

0.0 

875.0 

-778.9 

0.0 

*84 

1500,0 

-1294.4 

0.0 

800.0 

-855.8 

0.0 

800.0 

-850.9 

0.0 


a. Eigenvalues have no order relationship to value. 



















The slope of the pump curve may now be evaluated as 


ar _ d 33 l d 
A 1 B 1 S 0l 


where AjBj is a constant and Sqj is the pump speed. The compliance of the duct is 
simply 




(5-8) 


The resistance of the duct is 



-d 44 l d 

2W d 


(5-9) 


The inertance of the small fluid segment between the accumulator and the orifice is 
calculated as 



1 

°65 


(5-10) 


and its equivalent resistance is 



” D 66 L N 

2W n 


The accumulator parameters are determined similarly with 


C A = 
l a = 


_ 1 _ 

^78 

1 

®85 


(5-11) 


(5-12) 


(5-13) 
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and 


R A " ' D 88 L A • 


( 5 - 14 ) 


The Nonunique Equilibrium of the Sixth-Order Case 


The CTL-V test configuration has some unusual properties if the system is tested 
without the accumulator. The system no longer possesses a unique equilibrium, but is in 
equilibrium everywhere that the flows become equal. The system is stable with eigen- 
values, as reported earlier, but the system has infinite equilibrium conditions. The sixth 
order system is described as 



( 5 - 15 ) 


Now clearly 


AP 0S = (cjjAF FL - - AF 0S 


B 


( 5 - 16 ) 
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is zero if the two flows become equal, and 



dw fl af fl 


(5-17) 


determines the steady-state value of APqs for the nonzero AFpp. In the pump equation 


AP ODl = (^g) AFpL + D 2 AP ODl ' D l AF OS ‘ D 2 AP OI2 


(5-18) 


if APq D j is zero, then APqoi is related to the other variables as 

ap odi = ^ 2 ) ( D,AF ° S ’ D2AP ° 12 * (eg) AFpL 


(5-19) 


Proceeding, 


AF, 


OS 


= (i ) AI>0U1 "(~Ld DW ° S ) AF ° S (i) 


AP, 


012 


(5-20) 


which, if AFqs > s zero, similarly can be solved for APqqj as 

'2R D 


AP 


ODI 


“ L d((t-‘) 

\\ l D/ 


AP OI2 + I ' DW Os) AF 0S 


-D 


(5-21) 


The next equation is 


(=;) 


AP OI2 “ ( n~ K AF 0S * af OP2^ 


(5-22) 
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and finally 


AF 


0P2 


= (li) AP 0I2-(lJ DW0P2 ) AF ° P2 


(5-23) 


which yields a steady-state value for APqj 2 when all the flows are equal. Notice that by 
Eqs. (5-19) and (5-21) there are apparently two definitions of AP OD1 . Both are of the 
similar form 


AP ODl “ ap OI2 + C ' AF 


(5-24) 


where in Eq. (5-19) 


C' = 



(5-25) 


and in Eq. (5-21) 


e - 2 r d dw os 


(5-26) 


Interestingly enough for CTL-V at the operating point corresponding to the 1 00 percent 
engine power level, one finds that 



Therefore, any time that flows become equal, the system will be in equilibrium at that 
perturbed condition. This property is due to the particular values of the pump com- 
pliance, head rise characteristics, steady-state flow, and the duct resistance. 
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This result has no meaning in the context of the engine because the engine system 
has a closed fluid path around these elements, thus altering the overall system dynamic 
characteristics. Further, if the accumulator is added to the CTL-V system, then the sys- 
tem regains a unique stable equilibrium since the flow into the accumulator must go to 
zero in the steady state. This response poses no real problem to the technique of this 
research since the system is driven in an oscillatory fashion about the null, as is true of 
CTL-V itself. 


Results and Conclusions 

The most significant result is the ability of the observers or state reconstructors to 
follow the small signal nonlinear signal even though only an estimate of the system is 
known. Results in Tables 5-4, 5-5, and 5-6 demonstrate that even with estimates that are 
in error by large amounts, the robust observer provides reasonable estimates of the state 
for use in the estimation process. The nonlinear small signal values are the oscillations 
about the system operating point. The linear values are results from an analytic lineariza- 
tion of the nonlinear system equations. The first group of numbers is the linearized sys- 
tem Jacobian. This matrix is the analytic linearization of the nonlinear equations at the 
system operating point. The second group of numbers is the result of the estimation 
process at an instant of time, shown as the first number in the third group. The first two 
matrices may be compared by positions. The first row of eight numbers in the third 
group are the nonlinear states described in the first line separated by comas. The remain- 
ing rows are as described above. The constrained results of Table 5-6 refer to the method 
of parameter calculation. Parameter recalculations are permitted only for those elements 
that, due to model structure, are dependent. That is, accumulator parameters are not 
permitted to be a function of line or pump. The unconstrained estimation allows varia- 
tions as with any sensitivity technique. These differences may be observed in the tables as 
the estimate of the system Jacobian. These observers are not simply following the system 
but are reconvcrging after each reevaluation of parameters each 0.02 sec. Since the 
highest observer root is 200 Hz and the discharge has a resonance at 60 Hz, 2 to 4 sec is 
a very long run time. Either of these strategies work. The constrained method is similar 
to a steepest descent technique. The feedback gain must be small, for the constrained 
approach, to maintain computational stability. 

Figures 5-1 through 5-9 demonstrate the observer response. Notice that the pres- 
sure initial conditions are presumed known while there is error in the flow initial condi- 
tions. The observer response quickly eliminates the flow errors before the estimation 
process begins. This is a consideration in specifying a settling time. Figure 5-5 is a blowup 
of Figure 5-4 better displaying the observer response. The four second time histories of 
Figures 5-1 to 5-9 are to demonstrate that the observer response does not diverge over a 
long time interval, so that reasonable estimates of the unmeasured states are available for 
long time periods. 

Figures 5-10 and 5-1 1 show that the linear response is different from the small 
signal nonlinear response. Recall that the observer attempts to follow the small signal 
nonlinear response. 
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TABLE 5-5. OBSERVER PERFORMANCE AFTER 1 .36 sec (UNCONSTRAINED) 
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TABLE 5-6. OBSERVER PERFORMANCE AFTER 0.81 sec (CONSTRAINED) 
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Figures 5-12 through 5-17 are a sample of the estimation process response by 
matrix element. For completeness. Figure 5-14 shows the response of a zero element. All 
these results are for the eighth order, with accumulator case. 

Now examine results of the sixth order, without accumulator, configuration. 

Figures 5-18 through 5-23 show the response of the sixth order observer for the 
varying estimates of the parameters. Again the nonlinear small signal is being followed by 
a linear observer. Figures 5-24 and 5-25 show the response of the nonlinear small signal 
states contrasted with the linear response. The estimation response for the sixth order 
example is demonstrated by Figures 5-26 through 5-29. All of the responses were excited 

by a 10 Hz perturbation with an amplitude of 27.5 lb/in/. 

As was anticipated, the technique is fraught with sensitivity, numerical, and 
dynamic difficulties. The observer must be properly designed with respect to the observed 
system. Appropriate time intervals must be chosen to allow different dynamics to settle 
before the parameter estimation process begins. Sample rates must be chosen so that 
there is numerically sufficient change in variables, providing well conditioned matrices in 
the calculation process. Gains must be chosen properly to achieve adequate evaluation 
stability in the parameter estimation process. These considerations can be formidable, 
especially for systems possessing a large span in eigenvalues. Trial and error generally 
provide an adequate performance index to determine the times, gains, and sampling inter- 
vals. The technique possesses many shortcomings in terms of implementation and system 
dependence. However, the system has been demonstrated to work for the CTL-V con- . 
figuration, and the technique has many advantages in application to systems whose states 
cannot be directly measured. 
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Figure 5-15. Eighth order estimation 
response for D43. 


Figure 5-16. Eighth order estimation 
response for D33. 
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Figure 5-17. Eighth order estimation 
response for I>2\- 
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Figure 5-26. Sixth order estimation 
response for E> 2 J • 
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Figure 5-28. Sixth order estimation 
response for Dgg. 
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Figure 5-27. Sixth order estimation 
response for D 43 . 
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Figure 5-29. Sixth order estimation 
response for 
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APPENDIX A 


CTL'V SIMULATION WITH CROSS REFERENCE 
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034C 


7 

/S8E1E.85 > 

340.000 

n341 



Crt T 6 700 

‘ 341.000 

f)342 

752 


«Jo = 0 

342.000 

0343 



if usw.NE.uee to 7588 

343.000 

<1344 



DA<1/1)=C. 

344.000 

0345 



ca 7581 laO/KK 

345.000 

0346 

7581 


Oa( I / 1 ) =c * 

346.000 

0347 



DB 7582 l=4iNK 

347.000 

n346 

7582 


Da ( [ i c ) sC • 

348.000 

0342 



Da(1/3)«c» 

' 349.000 

ojSC 



Da (2j 3 ) sC» 

350.000 

n3Sl 



Cfj 7bfc3 laS/hN 

351 .000 

0352 

7583 


Da (1*35=q* 

352. OQO 

03S3 



CA { ?* 4 i =C* 

353.000 

o3 5 »< 



Dfi 7584 lafi.KN 

354.000 

03S5 

7584 


DA ( 1/4 ) aC* 

355.000 

o356 



Da ( 1*5 ) =o» 

356.000 



n357 


Ca(2/E)sG. 

157,000 

n35S 


Ca ( 5/ £ ) iC< 

.958.000 

n35S 


0A(7<5)aC* 

955*000 

036C 


CO 7585 I.IjA 

960*000 

036 l 

7585 

0ACl/6>=C* 

361 *000 

036£ 


CaC7/6)=C. 

362.000 

0363 


Oa ( a/ £ ) • C . 

363*000 

n36^ 


Co 7586 

964*000 

n365 

7586 

C a ! I / 7 ) » c • 

965*000 

0366 


1)0 7587 Uu4 

966*000 

0367 

75 37 

DA( 1 . 8 ) b C • 

367*000 

036S 


D A ( 6/3 ) »C* 

968*000 

0J6S 

738a 

DO 77C I s 1 / 1\ N 

369*000 

037C 


Co 76c t * 1 j N'N 

970*000 

0371 

760 

D ( I / j ) =>C ( T.j) +CA ( I/.. )*AON 

971 .000 

o37£ 

77u 

Continue 

372*000 

0373 


Go T6 758 

373.000 

0374 

790 

continue 

374.000 

n37fc 


IF (T.GT.TLM) call exit 

375,000 

0376 


Gfi t e 2oc 

976.000 

0377 

4 

FhR|>aT (/1F1G«3/(4E16.8) ) 

377.000 

037£ 

C ZEfifi 

I'Ef-aKV CALL 

378.000 

0379 

800 

CALL ZAF(KjM) 

.979.000 

038C 


go re lo 

.980.000 

0381 


EnO 

381.000 


gS 

v-d S 

l! 


^*3 



<1 

vo 





eRBSS 

REFEREaCF 
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1 



A 

13.0CC 

64 *0C0 

66*080 

67.000 

68*000 

69.(100 

70*000 

71*000 

72.000 


73*000 

74*0CC 

75*000 

76.000 

77 «0CQ 

78.000 

79*000 

80*000 

81.000 


82 • OGC 

83*008 

84*000 

85 .000 

87 * 8C0 

103.000 

135*000 

208*000 


AA . 

24*000 

48*008 

121*000 







ACN 

24. OCC 

48.0C 8 

371. OCO 







B 

13*CCC 

278*008 

286 * GCO 

P88.0C0 

289*000 

291 .000 

296*000 

?97 *000 


r 

13*0CC 

42*008 

48*000 

98.000 

198*000 

199*01)0 




rA ' 

42.0C0 

48*008 

82*000 

128. OCO 






rB 

42.00C 

48.008 

66 *CCO 

70. OCO 

122.000 





rD 

42 *0CC 

48*008 

77 *008 

126.000 






C8N 

42*008 

48 *ocr 

71 *0CO 

124 *0C0 






0 .. 

13*000 

42*008 

48*000 

1C3.000 

137*000 

138*00(1 

140*000 

208*000 

256*000 


220*000 

322*080 

324*000 

371 .000 






nA.. 

13*CCC 

304*088 

306*000 

322. OCO 

3?4 *000 

344.000 

346*000 

348*000 

349.000 


350*0CC 

352 *0C C 

353 *0C0 

355 *000 

356*000 

357 *000 

358*000 

359*000 

361.000 


362 *OGC 

263*000 

365*000 

367. OCO 

368*000 

371 .000 




OCl 

13.0CC 

122*008 

123*ocn 

124.000 

125*000 

126.000 

127*000 

128*000 

129.000' 

- . 

1 48 *OCC 

188*000 








80S 

13*008 

131 *008 

135 »0Cn 

142. OCO 

149*000 

189.000 




na3- 

13*000 

132 *088 

137*000 

138. OCC 

143*000 

150 -000 

190*000 



oD4 

13*000 

123 *0C8 

140*0cn 

144 .OCO 

151 *nco 

191 .000 




nE 

13*000 

254 *0C0 

256 • OCO 

259.000 






DT 

24 *OCC 

48 *008 

112*000 

114.000 

148. nco 

149.000 

150*000 

151*000 

164.000 

801- 

42*000 









nos 

42* OGO 









0 U 3 

42*000 









81 

13*000 

116*008 

122*000 

123 .000 

134*000 

125 .000 

126*000 

127*000 

’ 128.000 

. _ 

12S • OOC 

138*000 

154*000 

160. OOC 






82 

13*000 

117*008 

13b«CC0 

137.000 

155 < 000 

161 .000 




83. . 

13*000 

118*008 

137*000 

138 .000 

1F>6*000 

162.000 




84 

13*000 

119*008 

140*0C0 

157 .000 

163*000 





85 

114*000 

121. OCO 

158*000 

164 .000 

185. OCO 





F 

13. OOC 

25C *OCC 

256. CCO 







FO 

13*000 

252 * OCO 

259. COO 







FR 

13.000 

251* OCO 








frc 

13*000 

253*000 








fxh 

375*000 









p 

13*000 

313*000 

315*000 

324 .OCO 






T 

42*0C0 

48*0cn 

57*000 

58 .OCO 

62*000 

64 .000 

87*000 

94*000 

96.000 


. 98 * OCC 

103*008 

106.000 

ns. occ 

116. OCO 

117.000 

118*000 

119*000 

130.000 


131 *000 

132*008 

133*000 

135 .000 

137. nOO 

148 .000 

140 *0oo 

147*000 

148.000 


149*000 

15c*Oco 

151 .CCO 

154.000 

155*000 

156 .000 

157 *000 

160*000 

161.000 


162*000 

163*080 

179. COO 

180. OCO 

181*000 

182 .QOO 

183*000 

187*000 

188.000 


189*000 

1SC.0CO 

191 *CCO 

194.000 

195*000 

196.000 

198*000 

199*000 

202.000 


203*000 

2C8.0CC 

216* OCO 

218.000 

2?1 • OCO 

240.000 

241*000 

242*000 

249.000 


250*000 

251*000 

252. OCO 

253.000 

254*000 

256.000 

259*000 

260*000 

264.000 


265*000 

266*000 

270*000 

271. OCO 

£72. OCO 

274*000 

276.000 

278*000 

280.000 


281*000 

283.0C0 

284. OCO 

286*000 

£88.000 

289.000 

292*000 

293*000 

302.000 


304 * OOC 

306*000 

311*000 

313. OCO 

315*000 • 

320.000 

322*000 

324*000 

345.000 


346 * OCO 

347*000 

348. CCO 

351.000 

352*000 

354 .000 

355*000 

‘360*000 

361.000 


364*000 

365*000 

366 • OGn 

367.000 

369.000 

371.000 







CO 
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UN 

22 «QCC 

48. OCO 

136. OCO 

138. OCO 





1FS 

22 *00C 

48. ocr. 

2l2»ccn 

213.000 

Si 4*000 




TPS c 

22.000 

48. OCO 

320*000 

322*000 

334, noo 




IPS! 

22 » OOC 

4B*OCO 

iCS’OCn 

207. OCO 





r SVl 

22 »0C0 

4 8* OCO 

343.000 






“i 

42 »OCC 

48 * OCC 

63*000 

64.000 

87.000 

95.ono 

96*000 

98*000 


134.000 

135.000, 

137*000 

138*000 

l4o*nco 

197.000 

198*000 

199*000 


£48,000 

253. OCO 

260, OCO 

262.000 

£63,nco 

2&8.000 

263.000 

271.000 


£91 * OCO 

296*000 

237* ccn 

2SS.0CC 

3oo*noo 

303.000 

304*000 

306*000 


313*000 

315. OCO 

320«CC0 

322. OCO 

334 . nco 

370*000 

371*000 



S2.C0C 

246. OCf 

247. OCO 

548.000 

£63. noo 

285.000 

286*000 

388*000 


295 »OCC 

236. OCO 

342«CCO 






,'jK 

113. OCC 

145 • Qco 

146*000 

148.000 

l49.noo 

150*000 

151*000 

152*000 


154 *000 

155 • OCC 

156*000 

157.000 

l6Q.nco 

lftl.000 

162.000 

163*000 


ls>S»COC 

171* OCO 

172. CCo 

174. OCO 

175,000 

177.000 

178.000 


« 

12*CCC 

13.QC0 

42*000 

48*000 

ss*nco 

379*000 



KA 

24* OCC 

48. OCO 

60 *Qcn 

301*000 





KB 

6 C »CC5C 

. 237 >000 

238*000 

244*000 

3oi.nOQ 




kC 

12*000 

13. Ocr, 

96*000 

58*000 

13? *000 

138.000 



KK 

255 «00C 

256* OCC 

277*000 

278*000 

£87*nOQ 

288 . oon 

293*000 

894*000 


£97. OCC 

305. OCO 

3C6*CCo 

314.000 

315. nco 




kL 

24 .OCC 

48 .Occ 

235.000 






1 

97*0CC 

38*000 

284*000 

385. OCO 

287.000 

594,000 

295*000 


1 A 

12*000 

42*000 

48*000 

S3. OCO 

85* noo 

129*000 



I 0 

12*000 

42*000 

48*000 

71.000 

74.(500 

75,000 

124*000 

125*000 

i K 

24. OCO 

48*000 

2C6«0C0 






1 L 

3 2 'COO 

42. OCC 

48*000 

68 * 000 : 

69*000 

133.000 

142*000 

143*000 

1 H 

12*000 

42*000 

48*000 

80*000 

81*000 

127*000 



C 

282 * OCC 

283*000 

292*000 

379.000 





PE 

12*000 

13.000 

253*000 

306 . 000 

332*000 

334.000 



FK 

24 .000 

42*000 

48*000 

1S4.000 





vx 

12. OCC 

13. ocr. 

260*000 

278.000 

315.000 

322.000 

324*000 


NX!!' 

12* OCC 

13*000 

281*000 

231*000 

296. OCQ 

297.000 

306*000 

315*000 

kU 

56.CC0 

205. Oco 

2 C 6 . 0 CO 

2 C 8.000 

216.(500 

218.000 

281.000 

236*000 

Mv 

24 *000 

42*0co 

48*000 

57.000 

Io3.oCO 

115.000 

130*000 

134*000 


17S.QCC 

187. OCC ' 

137*000 

308*000 

£40*000 

247.000 

243*000 

255*000 


SbS-QGQ 

270*000 

276*000 

377.000 

£80*000 

282-000 

285*000 

287*000 


2S3.Q0C 

295. OCC 

302*000 

303*000 

3 o 5 .n 00 

31 1 .000 

312*000 

314*000 


347 *000 

351*000 

354*000 

369.000 

370* (ICO 




P 

13 * OCC 

27£*OCC 

274*000 

283*000 

£9l*nOQ 

296.000 

297*000 


pgpf 

42-000 

4 8 » OCO 

71*000 

1S4.0CO 





p0Pc 

42 « OCC 

48. OCO 







PFT 

42 *000 

48- Ocr 

123*000 






FT 

121 *000 

123*OCO 

142*000 

143,000 

144.(500 




POPS 

127*000 








RA 

42*000 

48. OCC 

85*000 

123. OCO 





rd 

42.000 

48. QCO 

72*000 

75.000 

lP4*noo 

135.000 



RL 

42.000 

48. OCC 

6 S* 0 co 

123.000 





RN 

42 .000 

48*0CC 

81-GC0 

127,000 





Rlf« 

121 .OCC 








Stfl 

42 .OCC 

48. OCO 

71 »ccn 

124,000 






103-000 

208.000 

281,000 

318.000 

289.000 

153.000 

168.000 


296.000 


144.000 


324*000 

147.000 

264.000 

231.000 

345.000 



T 

4e»0CC 

CROSS 
48 .Oco 

REFERekcf 

lc3*ccn 

tL^ 

£08 .OCC 

2 1 6 * 0 C o 

218. Con 

42 *CCC 

48 -OCO 

375. con 

U 

24 .COC 

48 .OCC 

121 »ccn 

X 

13.QCC 

42. OCC 

48.CCO 

xO 

13C*0CC 

2C3*0CC 

2C8*OCO 

13 .OCC 

188. OCC 

XCL 

13 >OCC 

189 .Oco 

253«0Cn 

XDM 

13'0CC 

191 .Oco 

242. CCO 

XDT 

13* OCC 

ISO »0CC 

242. ccn 

XL 

13 .OCC 

42. OCC 

46-CCn 

xLn 

1SS.OOC 

2C8.0CC 

£18. CCO 

13 * OCC 

58 «ocr 

1C3*CC0 

xt' 

13. OCC 

42«0CO 

46.CC0 . 

xS 

221 * OCC 

241 . OCC 

25c*Ccn 

13 » OCC 

42*0C0 

48.CCO 

XT 

13. COC 

42 «0CC 

48. CCO 

X2N 

195. OCC 

2C8.0CC 

21 2*000 

42. OCO 

48*0Cc 

6S » CCn 

X4K 

42*0CC 

48-000 

72-CCn 

XbN 

42 «0CC 

48 .OCC 

Sl-CCf! 

V 

13*0CC 

195‘Occ 

198-CCO 

VT 

13*0CC 

196. oco 

19S-CC0 

7AR 

71 

379* OCC 

1 3 * OOC 

148 »0CC 

154. ccn 

72 

13 *0C0 

149*000 

155* ccn 

73 

13*0CC 

15C*0CC 

156. ccn 

"74 

13*0CC 

151 *OCc 

157. CCO 

1 

id 

47.0CC* 

42 *0co 


21 .0CC4 

38C *Occ 


i i 

32C.CCC 

327 .OCO* 


1 lo 

S5.0CC 

S9-0C0* 


1 2 

322. OCC 

3£8. OCC* 


1 20 

94. OCC 

ICC. OCC* 


1 3 

324. QOC 

333 . OCO* 


1 9o 

102*0CC 

1C6.0CO* 


? 

25 • OCC* 

48 »0cr 


PCO 

108 »OCC* 

376. OCO 


3 

23.CCC* 

22-OCo 


30 

b 7 • OCC 

59 .Oco* 


300 

115. OCO* 

166.0CO 


310 

115 .OCC 

12C * OCO* 

lfc7. ocn 

315 

136.000 

138. OCC* 

320 

134. OCC 

14C-0C0* 


330 

13C.0CC 

141 .OCC* 


340 

147.00C* 

17C.0CC 

173* CCn 

350 

153* OOC 

1 6C • Qcc* 

3 60 

152.00C 

159. OCC 

165* CCO* 

370 

147. OOC 

166. OCO* 

380 

145 • CCC 

168* Oco* 


390 

lfafi.OOC 

171 «Ocn* 
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1C6.0C0 

H2. nco 

114. non 

221 .OCO 

375. nco 


58. OCO 

lc3.nco 

l u6. non 

216. OCO 

2i8.nco 

221.000 

252. OCO 

253. nco 

266. COO 

252. OCO 

266 . nco 


1C3.0C0 

106 «nco 

117.000 

221. OCC 

£51 • nco 

1C6.0CO 

2o3»nco 

208.000 

1C3.0C0 

1 i 5. nco 

157,000 

251.000 

265. nco 


58. OCO 

138. oCO 

2o3.oon 

1C3.0CG 

ic6.nco 

l i 8. con 

221 .OCO 

241 . nco 

250*000 

75.0C0 




3 


isa.ooo 

164*000 

185.000 

1 16 « OOC 

154*000 

160.000 

155.000 

'161.000 

181.000 

218.000 

163.000 

,221*000 

183*000 

208.000 

156*000 

2 6Q *CQQ 

162*000 

265*000 

182.000 


1 ^ 0*000 180'fico 
161 • oco m.nco 
162. OCO 182. nco 
163. OCO 183. nco 


176.000 



*t) S 

SI 

1® 
. a® 


3 Sl 

171. CCC 

qS 2 

174 .OoC 

953 

17S » 00 C 

364 

177. OCC 

395 

1«7 • OCC 

4 

377. OCC* 

4C 

O4-0CC. 

400 

197 * CCC 

4 1 0 

IS*. cr c 

415 

2c2«Cc c 

420 

207. OCC 

421 

212. OCC 

422 1 

213.0CC 

423 

214 .one 

424 

£15.000 

480 

288. CCC 

485 

c4C ►COC 

490 

23S.COC 

5 

a 8 » OCC* 

50 

62*000 

5C0 

2-.6.CCC 

510 

255 .OCC 

52C 

249 .OCC 

530 

264. OCC 

6 

216. OCC 

6 CC 

27C*0CC* 

610 

271. OCR 

62 0 

2 7C* OOC 

630 

273. OCC* 

63E 

£«i .one* 

640 

276* CCC 

650 

£ci3*Q(’G 

658 

2rtS . OCC * 

600 

296* CCC* 

670 

2 S 2 .CRC 

7 

lOfc.CCC 

700 

305. C«'C 

710 

3o3 *OCC 

720 

Sue » OCC 

755 

315 . OCC* 

756 

3 1 2 • C 0 C 

757 

3ll .CCC 

758 

315 .OCC* 

7581 

346 -CCC* 

758? 

318.0CC* 

7583 

352. OCC* 

7584 

3H5 • OCC* 

7585 

3al .OCC* 

7586 

265*000 

7587 

3l>7 *GCC* 

7588 

3l3 . CoR 


00 

w 


mtl9 FaR 10i '76 


fheSS RCf-ERErvCF ftF f'Alh 

171 *CCC* 

177 *Ccr* 

131 ‘OCR* 

137. OCR’* 

132 »OCR** 

6 3 -OCR 
2GC "OCR* 

2l 1 .OCR* 

20^ •OCR* 

pil.ocr* 

216 • OCR* 

218.0CR* 

£21 "OCR* 

2 17 * OCR 28000(1 223. OCC 234. nCO# 

24C " OCR* 

213‘Ocn* 

241 -OCR* 

27 ’ OCR 
65‘OCR* 

218« OCR* 

257. OCR* 

26 1 "OCR* 

267. OCR* 

224. OCR* 

3GC.0CR 
274 "OCR* 

273. OCR 275. OCR* 

277 .OCR 
i8C ‘OCR 
279* OCR* 

285 *0C0 B 2 6 • C C 0 230.0CC* 

237 .OCR 
295‘Occ 
258 »OCR* 

218 »OCR 226 • CCR* 

307 .OCR* 

3C8.0CR* 

309 .OCR* 

314. OCR, 

316. OCR* 

317 .OCR* 

373 .OCR 
345. OCR 
347. OCR 
351. OCR 
35l»0CR 
36C*QCR 
364 *OCR 
966 .OCR 
369 "OCR* 



00 

4 ^ 


iml9 I^aR 10* f 76 


5 


cRess 


759 

3 IS *OCC 

342 »0C0* 

760 

371*000* 

37C*ocr 

770 

369*000 

372*000* 

79q 

235«0CG 

236*000 

R 

221*000 

23C*0CC* 

800 

2C ■ OCC 

379*000* 

9 

24 *QCC 

54*000* 

99 

9S *OCC# 

97 *0C0 


REFEREnCF ftF pain 


244*000 268. OCO 


341*000 


374.000# 
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SIXTH AND EIGHTH ORDER DATA SETS FOR CTL-V EXAMPLE 
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- ? 2-1 _.o ... 1- ■ 

6 3 10 1 1C 

•00 1 » cl *5 ■> 6?*&j • C-GG5 * 
•0C 7 * *3.* .CCr.M6*5CC ** 

* 00 C iCl j - 1 • 4 * • c 1-3* • G C'l 1 H * 
635* * 635 * * *cCe# • C C 3 * 

• OGCl U4-.635.* .C-36* *COl> 

•?* ^ • * 3cO 

1 33* * 633*1.* 4 3 o • « j 633* 1* 

- M • 6fl j 633 • 1 # C • * (' •-* 

0* * • r * G • * *3* 

0 • * *3*0»jG»* 

0«*C*#C**G«* 

1 2 G ^ 1 9 S n * * 0 • / ( 'j • * 

147 7 • * 0 • ^ G * / C * * 

0 • j 

C • j ^ * Q * t ( • » 

G«*1»*C«*U»* 

• ...0 *>^#iO**C** - 

1 • * .C • t o * » c • > 

133* »A33*lj» it Su«^>j633 *1* 
•68>633*1* 1 H*S8 j*CC£* 

0 • f 

- 143-? **1 .06*143* jO«* 

CuC.iOiiQij 

C • j C_» j”32»8j 77 * .* 

0«j C»,/0.*0.* 
•143**0**“95*3j»111*4ji 
130 • » C * * 0 • * C * * 
C«-iG**32.*8*- 7.7_* j* . 

0. j 3.q3* jO»* ICOO* j 

0 » * C • J.Q aiOii 

-12S*i«^4«C3iO*jC** 

0 I J C «.* 0 • * C • j. . .. 

0 • j C • * 0 • * *9SC • i 
C*iC ( iOwO«j 

-125.jO**17*5i-2CC«^ 
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